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We study the graviton Regge trajectory in Holographic QCD as a model for high energy scattering 
processes dominated by soft pomeron exchange. This is done by considering spin J fields from the 
closed string sector that are dual to glueball states of even spin and parity. In particular, we 
construct a model that governs the analytic continuation of the spin J field equation to the region 
of real J < 2, which includes the scattering domain of negative Maldelstam variable t. The model 
leads to approximately linear Regge trajectories and is compatible with the measured values of 1.08 
for the intercept and 0.25 GeV“^ for the slope of the soft pomeron. The intercept of the secondary 
pomeron trajectory is in the same region of the subleading trajectories, made of mesons, proposed 
by Donnachie and Landshoff, and should therefore be taken into account. 


I. INTRODUCTION 


II. HOLOGRAPHIC QCD MODEL 


The Pomeron plays a crucial role in QCD Regge kine¬ 
matics, for processes dominated by exchange of the vac¬ 
uum quantum numbers. This includes elastic scattering 
of soft states at high energies and low momentum trans¬ 
fer. The corresponding amplitude exhibits a universal 
behavior explained within Regge theory [T], 

M(s,t) «a(t) = 1.08-h 0.251, (1) 


We will consider the Holographic QCD model proposed 
in the works [HHH based on gravity plus a dilaton field. 
We shall be working in the string frame because the 
Regge trajectory we are interested in is made of funda¬ 
mental closed string states. As usual, the scalar field 
$ = d)( 2 ;) and the dual geometry has metric 

= gabdx’^dx^ = [dz^ -\- r]apdx°'dx^) , ( 2 ) 


in GeV units and for some function j3{t) that depends on 
the scattered states. A precise computation of the values 
of the intercept (ao = 1.08) and slope {a' = 0.25 GeV“^) 
is beyond our current understanding of QCD, since long- 
range strong interaction effects are important. 

The gauge-gravity duality is a new tool to unveil QCD 
strongly coupled physics [2]. In particular, the Pomeron 
is conjectured to be the graviton Regge trajectory of the 
dual string theory [3]. This fact has been explored in 
diffractive processes dominated by Pomeron exchange, 
like low-x deep inelastic scattering (DIS) [IHS], deeply 
virtual Compton scattering [7] , vector meson production 
[H] and double diffractive Higgs production [3]. 

Consider for instance the case of low-cc DIS. One ob¬ 
serves a rise of the intercept jo from 1.1 to 1.4 as Q grows, 
where Q is the momentum scale of the photon probe. The 
conventional approach is to start from the perturbative 
BFKL hard pomeron [10], which still exhibits conformal 
symmetry. Introducing a cut-off, one explains the ob¬ 
served rise of the structures functions and even the run¬ 
ning of the intercept, provided the cut > 4 GeV^ is 
imposed in the kinematics m- However, dual models 
that also start from a conformal limit and introduce a 
hard wall cut off in AdS space, give even better hts to 
data, without imposing any restriction in the kinematics 
[6|. This is a strong motivation in favor of treating soft 
pomeron physics using the gauge/gravity duality. 

This letter builds a soft-pomeron phenomenology in 
Holographic QCD. More concretely, we show the Regge 
theory for spin J exchanges in the dual geometry leads to 
the behavior 0 for the amplitude between soft probes. 


where r]a 0 is the Minkowski boundary metric. In the 
string frame the corresponding action is 



d^xy/—ge 


R + A{d<^f + V 


(3) 


with 


U = e 
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3 ) 


(4) 


The field 4> is the dilaton without the zero mode that 
is absorbed in the gravitational coupling k. The field 
equations arising from ([^ take the form 

1 dV 

Rab + 2VaV6$ — -7^ gab =0, 

4 

2v 24.-4(V$)^-bU-f =0. (5) 

4 dd) 

The superpotential W ($) is fixed phenomenologically by 
demanding that the model reproduces basic QCD data, 
such as beta function, heavy quark/anti-quark linear po¬ 
tential and glueball spectrum. In this work we take the 
Background I of miisi where 



:boX 


1-b 


{2bl -b 36i) log(I -b A^) 1 ^ 
I8a 


( 6 ) 

A = e'^ and the length scale L fixes the units. 

The’t Hooft coupling of the dual Yang-Mills theory A 
is fixed by A up to a multiplicative constant, i.e. A = cqA. 
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FIG. 1: Running coupling ots vs. energy scale. The red point 
is a<,(1.2GeV) = 0.34. 

For the model considered in this work the constants in 
([^ are given by 


III. POMERON IN HOLOGRAPHIC QCD 
A. Graviton 

Since we are interested in the graviton Regge trajec¬ 
tory let us start by considering perturbations to the back¬ 
ground in the string frame. We shall write the metric and 
dilaton, respectively, as 

ffab + hab , -b . (10) 

It is then a mechanical computation to obtain from (§ 
the linearised equations of motion for the perturbations 
hab and tp, 

'^^hab — + '^a'^bh -b 

-bdV^V(ad> hb)c + 2V'^$ (2V (^ahb)c ~ V chab) 

-4VaVb(^+igabR"($)¥> = 0, (11) 

vV+^R'(<i>)<^ + 5p"($)(^ 

-4V^ • {2\7^hab - ^ah) 

-/l“''VaVb$ -b 2/l“''Va$Vb$ = 0 , (12) 


foo = 4.2, 


bi 

% 


51 _ 3 

m ’ “ “ 16 ■ 


(7) 


The model has an additional integration constant that 
can be related to Aqcd, via the identification of the en¬ 
ergy scale and warp factor, \ogE = A{z) — |4’(2:). As 
shown in min], the UV behavior of the superpotential 
^ leads to a beta-function 


/3 


dX 

d\ogE 


= -boX^ - biX^ + ... . 


( 8 ) 


This is consistent with the two-loop perturbative beta 
function in large-N Yang-Mills 


^ = -boX^ - biX^ , 


- _ 2 11 5i _ 51 
° “ 3 (47r)2 ’ “ m ■ 


(9) 


if we take cq = bo/bo. This fixes the second parameter 
in 0. The others parameters are fixed by the IR con¬ 
straints coming from confinement and asymptotic linear 
glueball spectrum and lattice QCD. 

Given that all of the parameters are already fixed at 
this point, one may ask how the field theory coupling runs 
with energy. Setting Ac = 3 the QCD running coupling 
can be identified with as = A/(127r). Figure[^shows how 
Us runs with the energy scale in the model, giving 0.34 
for the value for E = 1.2 GeV, which is very close to the 
experimental value 0.35. 

We can recover the conformal limit by considering the 
parameters 5o = &i = 0. Then we can set <f> = 0, 
the superpotential becomes the cosmological constant 
— 12/L^, and the metric becomes that of AdS space, i.e. 
A{z) = \n{L/z). 


where the covariant derivatives and Riemann tensor refer 
to the background and h = h^. Field perturbations will 
be classified according to the S'0(l,3) global symmetry 
of the background. Thus we shall decompose the metric 
perturbations hab as 

hal3 = "b “b (49 q,i9^ — Jl + Tjaph , 

hzz , hza — -b daS • (^8) 

As usual transverse and traceless (TT) tensor fluctua¬ 
tions, transverse (T) vector fluctuations and scalar fluc¬ 
tuations decouple. Moreover, since we are interested only 
in the TT metric fluctuations, we do not need to worry 
about mixing of perturbations. It is then simple to see 
that 0 gives for the TT metric fluctuations 

-b = 0. (14) 

The term with the dilaton arises from the usual coupling 
—28'^^'^chab for metric fluctuations in the string frame; 
the other term comes from the coupling to the Riemann 
tensor Racbah^^"^,- T]apr]f,^) 
and Razfiz = —Ae’^^Vad- In the case of pure AdS space, 
A{z) = ln(L/z), so ( [l4| simplifies to 

{^^-m^)hlj = 0, (15) 

with {LmY = —2, as expected for the AdS graviton. 

B. Dual spin J field 

We shall consider the exchange of twist 2 operators of 
Lorentz spin J formed from the gluon field [23] 

Oj ~ tr , 


(16) 
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FIG. 2: Expected form of the A = A(J) curve (in blue). 


where D is the QCD covariant derivative. The dimension 
of the operator Oj can be written as A = 2 + J + 7 j, 
where 7,/ is the anomalous dimension. In free theory the 
operator has critical dimension A = 2 + J. 

Knowledge of the curve A = A( J) is important when 
summing over spin J exchanges, since this sum is done 
by analytic continuation in the J-plane, and then by con¬ 
sidering the region of real J < 2. Figure summarizes 
a few important facts about the curve A = A(J). Let 
us define the variable 1/ by A = 2 -|- w, and consider the 
inverse function J = J{v). The hgure shows the pertur¬ 
bative BFKL result for J(;z), which is an even function 
of v and has poles ai iv = \. Beyond perturbation the¬ 
ory, the curve must pass through the energy-momentum 
tensor protected point at J = 2 and A = 4. We shall 
use a quadratic approximation to this curve that passes 
through this protected point, 

J{v) « Jo - 4X> = 2 - Jo . (17) 


The use of a quadratic form for the function J{v) is 
known as the diffusion limit and it is used both in BFKL 
physics and in dual models that consider the AdS gravi¬ 
ton Regge trajectory (see for instance H)- 

Consider now the spin J field dual to the twist 2 op¬ 
erators (1^. For pure AdS this field obeys the equation 

(V^ - m2) = 0 , {Lmf = A(A - 4) - J, (18) 


where L is the AdS length scale. Note that this field is 
symmetric, traceless and transverse iV^hba^ -aj — 0). 

To consider the spin J field in a general background 
of the form ([^, we need again to do a decomposition 
in S'0(l,3) irreps. The propagating degrees are de¬ 
scribed by components ha^-.-aj, since the other compo¬ 
nents hz-.-zoLf-aj (* > 2) are fixed by the transversality 
condition. Thus we need to dehne the equation of motion 
for ha^...aj. Of course we do not know its form for the 
dual of QCD, but follow a phenomenological approach. 
We shall require that such equation is compatible with 
the spin 2 case (141, since in that case it must reduce 


to that of the graviton, whose dual operator has pro¬ 
tected dimension. Moreover, we require the coupling to 
the dilaton to be that of closed strings in the graviton 
Regge trajectory arising from the term 
Finally, we require the equation to reduce to (181 in the 


conformal limit (constant dilaton). This leads to the fol¬ 
lowing proposal 


V2-2e"2^$V,,- 




( 19 ) 

where here L is a length scale parameter. It is trivial 
to verify that setting J = 2 (and A = 4) this equation 
reduces to the the graviton equation (14). Similarly, set¬ 


ting A{z) = ln(L/z) and $ = 0 we recover the spin J AdS 
equation (18). The dilaton term arises from considering 


tree level closed strings, which is justified since we work 
at large N. We expect that there will be more terms in 
this equation arising from other curvature couplings and 
derivatives of the dilaton held. Assuming the equation is 
analytic in J, these terms should be proportional to J—2, 
so that they are absent for J = 2. Notice that there can 
be such terms still at the level of two derivatives, that is 
terms proportional to 


„-2A 




e-2^$2. 




( 20 ) 


which also vanish in the conformal limit. Terms with 
higher derivatives will appear in a a' j!?• expansion. As 
already stated, we shall follow a phenomenological ap¬ 
proach and use the simple form (19) to describe the fluc¬ 


tuations of the spin J held in holographic QCD. 

We will be interested in the continuation of (19) to the 


unphysical region of J < 2. It is here that we will use 
the diffusion limit (fTrl), writing in (IT^ 


( 21 ) 




with Is a length scale set by the QCD string. Notice that 
we are hxing Ig to a constant determined by IR physics, 
but in fact it should depend on energy scale, since the 
curve A = A( J) in Figure should vary with energy 
scale, keeping its general shape. However, for the soft- 
pomeron this should not matter [ 21 ]. We leave Ig as a 
phenomenological parameter to be hxed by data. 

In the Regge limit we are actually interested in the 
-I- • • • -I- component of (19). To hnd the solution write 


h+...+ {z,x) = , 


( 22 ) 


where q ■ x = riap<f‘x^ and we set g_ = 0 in the Regge 
limit. Then, a computation shows that (19) reduces to 
the Schrodinger problem 


(23) 


-f [/(z) ) il){z) = til){z ), 


U{z) = ^- + , (24) 

4 

with t = —q2. The energy spectrum for each J quantises 
t = tn{J), therefore yielding the glueball masses. 
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C. t-channel spin J exchange 


Next consider the elastic scattering of QCD hadronic 
states of masses mi and m 2 . We write the incoming 
momenta ki, k2 and the outgoing momenta ks, ki in 
light-cone coordinates (- 1 -,—,_L) as 


ki=(ys,^,0^, ks 
fc2 =(^,V^,0), ki 


fml + ql !- \ 


where we consider the Regge limit s'^ t = —<4- 

Each hadron is described by a normalizable mode 
Ti{z,x) = where ^3 = v\ and Vi = uj. The 

hadrons we consider are made of open strings. Then the 
coupling of each hadronic field to the spin J closed string 
fields has the form 

nj j d^x^e-^K,...a,T^‘^^ ■ ■ - . (26) 


Notice that in principle different types of hadrons will 
have a different coupling kj. The transverse condition on 
the spin J field guarantees that this coupling is unique up 
to derivatives of the dilaton field, which are subleading 
in the Regge limit. 

The amplitude for mi m 2 —t mi m 2 scattering through 
exchange of a spin J field in the t-channel may now be 
computed in the dual theory. In the Regge limit we have 



0 1 2 3 4 z 


FIG. 3: Effective potential for different values of spin J. The 
first 2'^'^ glueball states are also shown. 


is the Fourier transform of 
Gj{z,z',l±) 

- J d'w'^d'w~Ili __ {z,z',w), (31) 


where w = x — x' = {w'^,w~,l±) and l± = x± — x'j_. 

From the H-h,-component of (28 1 , as defined 

by (191, it follows that Gj{z, z',l±) is an Euclidean scalar 
propagator in the three-dimensional transverse space of 
the dual scattering process (da;+ = dx~ = 0 in ([ 2 )), i.e. 


Ajih) = -kjk'j j d^Xd^X'^44fe-^-^' 

(Tiai^Ts) n-’+ '+(X,X')(T^9(/T^) , (27) 

where X = {z,x) and X' = (z',x') are bulk points and 
fields with a prime are evaluated at X', e.g. = $(z'). 
We use this notation throughout. We expect the spin J 
field propagator to obey an equation of the type 


□3 - ^ _ 2 A 4 j 

A(A - 4)] ^ -e^^SAx, x') , (32) 


where here x = {z,x±) and x' = (z',x(j_). Writing 




Gj{z,z\t) = e 


A(z) 


i’iz), 


(33) 


{vnu...ajM-bjix,x') = 

9ai{bi ''' 9\aj\ bj)^ 5 (Al, X ^ traces , 


the homogeneous solution to ( |32[ ) is exactly given by the 

~24|). Moreover, using 


Schrodinger problem of (231 and 
(28) J2ni’n{z)'ilj4iz') = Hz —zH, we conclude that 


for some differential operator T). We are interested in the 
-h • • • -h, — • • • — component of this equation, for which the 
differential operator T) can be read from (19). 

Some algebra shows the amplitude (27) simplifies to 

-4j(s,f) s J dzdz'< 


^/g3^+3T'-$-3>' 


l^^ll 


J -1 


^ Gj(z,z',t), (29) 

where V is the boundary volume. The function 

Gj{z,z',t) = j dH^e-^^^-^^Gj{z,z',li_), (30) 




(34) 


Note the eigenvalues and functions ipn depend on J. 


D. Regge theory 


We will sum all even spin J exchanges with J > 2 
using a Sommerfeld-Watson transform 


f E G+(-T) 

.7>2 


TT f dJ + { — sY 
2 J 2TTi sin( 7 rJ) ’ 


( 35 ) 
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which requires the analytic continuation of the amplitude 
Aj{s, t) to the complex J-plane. Then, the amplitude for 
the exchange off all even spin J fields becomes 

A{sA) = iV f dzdz'e^^^~''^'^vi\‘^\v2\‘^^Xn, (36) 

n 

where x„ = Xn{z, z', s, t) is given by 

TT f dJ S'’ + {—sY Kjk'j 

2 J 2TTi sin(7rJ) 2-’ 

^-(J-1)(A+A')MY'4’Yz') 


We assume the J-plane integral can be deformed from the 
poles at even values of J, to the poles J = jn{i) defined 
by tn( J) = t- In the scattering domain of negative t these 
poles are along the real axis for J <2. Thus we can write 


Xn = 


2 V 2 




2i" 


.-Ur.-k)A+A')dj 




(38) 


where (t) and we remark that the wave functions 

ipn are computed at J = jn{t)- It is clear that for large s 
the amplitude ( [M| ) will be dominated by the Regge pole 
with highest j„ in accord with the Regge behavior Q . 

We now specify to the model considered in this letter, 
which is determined by the effective Schrodinger poten¬ 
tial ( [^ . Since we are interested in the region J < 2, 
we can use the model introduced in (21) for the curve 
A = A( J). Figure|^shows the potential for several value 
of J. The energy levels for J = 2 are shown and compute 
the mass of the spin 2 glueball masses. As J decreases 
the energy levels will eventually cross the zero energy 
value. This will be the value of the intercept for the n-th 
Reggeon. Figure]^ shows the curves jnit), which clearly 
show that n = 1 is the leading Regge pole. The curves 
are approximately straight so we can also define a Regge 
slope. Note that the model considered in this paper al¬ 
lows us to investigate the region of real J < 2 and find 
the Regge poles. This differs from previous approaches 
based on Regge trajectories (see for instance m)- 


IV. RESULTS 

Finally we can test to which degree we are reproduc¬ 
ing QCD physics. We consider first the leading Regge 
pole. We vary Ig, introduced in (pTj ), to fix the Pomeron 
intercept to the value given in [1], as an optimal fit for 
total cross sections. For the value 4 = 0.178 GeV“^, 
and independently of our choice of Aqcd, we obtained 
ao = 1.08. The value of the slope is then fixed by 
the choice of Aqcd- We obtained ct'AgQjy = 0.018. 
If we fix Aqcd = 0.292 GeV as in [TH [T3], such that 
the first glueball mass TOo++ = 1.475 GeV, one obtains 



FIG. 4: The first Regge trajectories that result from solving 
the Schrodinger problem for discrete values of J. 


a' = 0.21 GeV If, on the other hand, we require 
the measured value of a' = 0.25 GeV we ob¬ 

tain Aqcd = 0.265. This is consistent with having the 
2++ glueball of the Pomeron trajectory with a mass of 
1.9 GeV, which is a known possibility [13 [HI- 
Let us remark that we could fix Ig to reproduce the 
intercept obtained in lattice simulations of SU(?>) pure 
Yang-Mills [T5]. In this case, for Ig — 0.192 GeV“^ one 
has Qfo = 0.93. Then, setting A = 0.292, which is fixed 
to reproduce mQ++ = 1.475 GeV of the same lattice sim¬ 
ulations, we obtained a slope a' = 0.25 GeV“^. This is 
exactly the slope obtained by the lattice simulations m- 
For the second pole we obtained an intercept of 0.433, 
which is consistent with the value used in [T]. We ran 
fits to pp total cross section data [20] and found that 
the second pole is necessary and needs to be in a narrow 
range of « 0.35 — 0.55. We determined this range by 
fitting an expression of the form 

a = go{a'3^° + gi{a'sYA (39) 

using go and gi as parameters, and varying ai. We fit 
this to pp scattering data with > 10 GeV. The above 
range is fixed by the requirement that Xd o f order 1 
or less. Our results can be seen in figurej^ As can be seen 
there using just the leading Pomeron exchange fails to fit 
the data satisfactorily. The second pole in [T] corresponds 
to several degenerate meson trajectories, while here it 
represents a next-to-leading glueball trajectory. Thus, 
our work points to the possibility that in this range there 
is a glueball trajectory as well. In fact, at least some of 
the /2 states are known to correspond to glueballs (see 
[21] and references therein for recent results). 


V. CONCLUSION 

Soft pomeron physics is still beyond the current an¬ 
alytic understanding of QCD. The best one can do at 
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Total Cross Section 



FIG. 5: A fit to pp total cross section data using the exchange 
of the first two Regge poles in our model. The green line 
represents the leading Pomeron exchange, and fails to fit the 
data at moderate values of 

weak coupling is to start from the BFKL approach and 
then introduce the running of the coupling, therefore 
breaking conformal symmetry. As a consequence, the 
branch cut of the BFKL pomeron becomes a set of poles 
in the J-plane m- This approach can be used to fit 
DIS data for hard scattering, keeping a very large num¬ 
ber of poles. However, it is not applicable to the case of 
soft probes. In general we expect to have a description 
of soft pomeron exchange as a Regge pole, in agreement 
with the phenomenological approach pioneered by Don- 
nachie and Landshoff [1] . Such description was proposed 
in [3] , based on scattering of closed strings in a dual con¬ 
fining background. In particular that work anticipated 
that for confining theories with a negative /3 function the 
pomeron, described as the graviton Regge trajectory, be¬ 
comes a Regge pole. Our work confirms this expectation 
by extending the holographic QCD model of m to scat¬ 
tering processes dominated by soft pomeron exchange, 
bringing a new insight to soft pomeron physics. 


Let us finish with a caveat and two open questions. It 
has been claimed that a soft Pomeron pole is not enough 
to describe the new LHC data [52]. This is somewhat 
expected, since it is known that at very high energies such 
a Regge pole would violate the Froissart-Martin bound, 
and other effects need to be included, for example multi- 
Pomeron exchange. However, this does not invalidate the 
great experimental successes of soft Pomeron exchange 
up to LHC energies, as well as the necessity to understand 
the subleading trajectories. 

The first question concerns the relation between hard 
and soft pomerons. Recent studies in gauge/gravity du¬ 
ality reproduce a plethora of low-a:: processes using the 
graviton Regge trajectory as the dual trajectory of the 
QCD Pomeron mai. In these cases one observes a run¬ 
ning of the intercept with the size of the probes. It 
would be very interesting if we could embed these re¬ 
sults within the present model, therefore unifying both 
pomerons. Another question is related to the spectrum 
of the spin J field, at integer values. It would be very 
nice to reconstruct the spin J equation in this domain, 
such that it reproduces perturbative QCD results. 
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In the singlet sector there is also the twist 2 quark opera¬ 


tor Oj ~ %i)TaiDo 


■ Daj’ip, however we are considering 


processes dominated by exchange of the gluon field. 

For instance, the approximation (211 misses the dimen¬ 
sions of the operators with J > 2 in the free theory limit. 



